Rapid Cooling of the Neutron Star in Cassiopeia A 
Triggered by Neutron Superfluidity in Dense Matter 
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We propose that the observed cooling of the neutron star in Cassiopeia A is due to enhanced 
neutrino emission from the recent onset of the breaking and formation of neutron Cooper pairs in 
the 3 P2 channel. We find that the critical temperature for this superfluid transition is ~ 0.5 x 10 9 K. 
The observed rapidity of the cooling implies that protons were already in a superconducting state 
with a larger critical temperature. Our prediction that this cooling will continue for several decades 
at the present rate can be tested by continuous monitoring of this neutron star. 

PACS numbers: 97.60.Jd, 95.30.Cq„ 26.60.-c 



The neutron star in Cassiopeia A ("Cas A" hereafter), 
discovered in 1999 in the Chandra first light observa- 
tion [l| targeting the supernova remnant, is the youngest 
known in the Milky Way. An association with the histor- 
ical supernova SN 1680 [2j gives Cas A an age of 330 yrs, 
in agreement with its kinematic age Q. The distance 
to the remnant is estimated to be 3.4lg' 3 kpc The 
thermal soft X-ray spectrum of Cas A is well fit by a 
non-magnetized carbon atmosphere model, with a sur- 
face temperature of 2 x 10 6 K and an emitting radius of 
8—17 km [jj. These results raise Cas A to the rank of 
the very few isolated neutron stars with a well determined 
age and a reliable surface temperature, thus allowing for 
detailed modeling of its thermal evolution and the deter- 
mination of its interior properties @ . 

Analyzing archival data from 2000 to 2009, Heinke & 
Ho 0] recently reported that Cas A's surface tempera- 
ture has rapidly decreased from 2.12 x 10 6 to 2.04 x 10 6 K 
[H, [j| . This rate of cooling is significantly larger than ex- 
pected from the modified Urea ( "MU" ) process or 
a medium modified Urea [l2j . It is also unlikely to be due 
to any of the fast neutrino (V) emission processes (such 
as direct Urea processes from nucleons or hyperons, or 
v — emission from Bose condensates or gapless quark mat- 
ter) since the visible effects of those become apparent over 
the thermal relaxation timescale of the crust [l3j], i.e., 
30 — 100 yrs, much earlier than the age of Cas A, and ex- 
hibit a slow evolution at later times. We interpret Cas A's 
cooling within the "Minimal Cooling" paradigm [lj[ and 
suggest it is due to the recent triggering of enhanced neu- 
trino emission resulting from the neutron 3 P2 pairing in 
the star's core. Our numerical calculations and analytical 
analysis imply a critical temperature Tc — 0.5 x 10 9 K 
for the triplet neutron superfluidity. 

The essence of the minimal cooling paradigm is the a 
priori exclusion of all fast ^-emission mechanisms, thus 
restricting ^-emission to the "standard" MU process and 
nucleon bremsstrahlung processes llj. However, effects 



of pairing, i.e., neutron superfluidity and/or proton su- 
perconductivity, are included. At temperatures just be- 
low the critical temperature T c of a pairing phase tran- 
sition the continuous breaking and formation of Cooper 
pairs [15j, referred to as the "PBF" process, results in an 
enhanced neutrino emission. Calculations of T c for neu- 
trons, T cn , in the 3 P2 channel relevant for neutron star 
cores, are uncertain due to unsettled interactions and 
medium effects which can either strongly suppress or in- 
crease the pairing [Tt| . Consequently, predictions range 
from vanishingly small to almost 10 11 K 14j |. The pair- 
ing gap is density (/?) dependent, and the resulting T cn (p) 
commonly exhibits a bell-shaped density profile. Assum- 
ing the neutron star has an isothermal core at tempera- 
ture T, the phase transition will start when T reaches, at 
some location in the star, the maximum value of T cn (p): 
Tc = max (p). At that stage, neutrons in a thick shell 
go through the phase transition and as T decreases, this 
shell splits into two shells which slowly drift toward the 
lower and higher density regions away from the maximum 
of the bell-shaped profile. If the neutron 3 P2 gap has the 
appropriate size, ^-emission from the PBF process is an 
order of magnitude more efficient than the MU process 
(see Fig. 20 of or Fig. 2 of 0]). 

Implications of the size of the neutron 3 P2 pairing 
gap were considered in [3] with the result that, for 
Tc < 10 9 K, a neutron star would go through the pairing 
phase transition at ages ranging from hundreds to tens 
of thousands of years, accompanied by a short phase of 
rapid cooling. This phenomenon, illustrated in Fig. 6 of 
[HI, closely resembles earlier results, e.g., Fig. 1 of [3] 
and Fig. 8 of 0. In Fig. HJ we show similar results of 
models with the likely surface chemical composition of 
Cas A from observations @, Q ■ The dotted curves show 
that rapid cooling is strongly dependent on Tc, a value 
~ 0.5 x 10 9 K being favored by the present age of Cas A. 

We now offer physical insight into the evolution of 
Cas A based on an analytical model developed in [llj . 
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FIG. 1. Red-shifted effective temperature vs. time, of an 
isolated neutron star for various values of the maximum neu- 
tron 3 P-2 pairing critical temperature, Tc (in units of 10 9 K). 
The circled star shows the value for Cas A. Results shown 
are from numerical solutions of the general relativistic energy 
balance and heat transport equations for a 1.4 Mr star built 
from the APR equation of state [2(|. Refs. 0, [3| describe 
the microphysical inputs used. These results are similar to 
those in Fig. 6 of [Tj|, but with the difference that the up- 
per layers of the envelope contain light elements of total mass 
AAfnght = 5 x 10 M©, compatible with the observation of 
carbon at the star's surface [1, 01- The rapid cooling at ages 
~ 30-100 yrs is due to the thermal relaxation of the crust. 
Thereafter, the stellar interior is isothermal. 



Assuming the star's core to be isothermal at a temper- 
ature T, a state reached within a few years after birth, 
and that the dominant cooling occurs through ^-emission 
with a luminosity L v , which is correct for up to several 
tens of thousands of years, energy balance gives 



r dT - t 

at 



(1) 



where Cy is the star's total specific heat. Neglecting the 
effects of pairing [2l| , one can write Cy — Cg Tg with 
Cg ~ 10 39 erg K" 1 and Tg = T/(10 9 K). Similarly, L v 
from the MU process can be written Lmu — Lg Tg with 
Lg ~ 10 40 erg s . Once triggered, neutrino emission 
from the PBF process is more efficient than the MU pro- 
cess; i.e., I/pBF = JLmu with / ~ 10 [22J. Values of Cg 
and Lg can be found in [ll[ or from Figs. 11 and 20 of 
[l3 | . Considering ages large enough such that T is much 
smaller than the initial temperature, Eq. ([T]) gives 



T=10 9 K(r MU A) 



1/6 



(2) 



with a MU cooling time-scale r MU = 10 9 Cg/QLg ~ lyr. 
The temperature T will reach Tc at a time to — 
t mv /Tqq, If this happened recently in the core of Cas A, 
i.e., at tc — (0.5 — 0.9) x 330 yrs, one deduces that 

T c = maxT cn (p) = 10 9 K(r M uAc)^ ~ 0.5 x TO 9 K. (3) 



FIG. 2. Schematic evolution of the internal temperature T of 
a neutron star (solid curves). Left panel: Initially, T follows 
the upper modified Urea (MU) trajectory, Eq. p]l. until T 
reaches Tc at time tc when the pair breaking and formation 
(PBF) process turns on. Thereafter, T rapidly transits, on 
an initial time scale r TR = tc/f, Eq. Q, toward the lower 
PBF trajectory, Eq. ((6|. Empirically, the rapid transition 
lasts up to a time < 2tc- Right panel: Evolution using a 
superconductivity suppressed modified Urea (SMU) rate. For 
the transit to start at a time t' c w tc the trajectory requires 
a T'c larger than Tc of the left panel. The early transit has a 
shorter time scale t^, r =t' c /f with /' > /, but the late time 
evolution is almost identical to that of the left panel. The left 
(right) panel corresponds to models of Fig. Q] (Figs. [3] and [4]). 



After that moment, L v rapidly increases, by a factor /, 
and the solution of Eq. |T]) switches to 



T 
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with a transit time scale r TR given in terms of the PBF 



time scale r D 
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With an / of the order of 10, one can naturally expect a 
transit time scale of a few decades in the case of Cas A. 
When t ^> tc, the solution in Eq. @ has the same form 
as that in Eq. ^ , but with the shorter time scale t pbf : 



T = 10 9 K(r PBF A) 1/6 = 10 9 K(r MU //i) 
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The above evolution is depicted schematically in the left 
panel of Fig. O 

For comparison with observations, we convert the in- 
ternal T to the observed effective T e using (see, e.g., |ll| ) 



T e ~ T e0 (T/10 8 K)' 3 K , 



(7) 



where T e0 - 10 6 K and (3 - 0.5. The evolution of T e 
is hence similar to that of T, and the internal cooling 
curves of Fig. [2] map onto analogous models of Figs. [TJ 
131 and |U The scale T e0 and the exponent (3 in Eq. (0) 
both depend on the chemical composition of the envelope. 
The presence of light elements, e.g., H, He, C, and/or O, 
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increases T e o and reduces j3 compared to the case of heavy 
elements, e.g., Fe. depending on the total mass AAfi; g ht 
of light elements [23| . 

Using Eq. (0, the slope s = d\og 10 T e /d\og 1Q t of the 
transit cooling curve from Eq. (Q| is 



d\og w T 
d\og w t 
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whereas the slopes of the asymptotic trajectories, 
Eqs. ([2j and ©, are both s = -(3/6 ~ -1/12. As long as 
t is only slightly larger than tc, the transit slope is larger 
than those of the asymptotic trajectories by a factor ~ /. 
The observed slope over a 10 yr interval is s b s — —1-4. 
Note, however, that the model "0.5" of Fig. [T]does not 
exhibit such a large slope. We are thus led to investigate 
the origin of the rapidity of Cas A's cooling. 

Several factors influence the rapidity of the transit 
phase. Firstly, Lpbf depends on the shape of the T cn (p) 
curve. A weak p dependence, i.e., a wide T cn (p) curve, 
results in a thicker PBF neutrino emitting shell and a 
larger Lpbf than a strong p dependence. Secondly, the 
T dependence of T e , i.e., the parameter fj in Eq. ([7]), 
also affects the slope in Eq. ([5]). Thirdly, protons in the 
core will likely exhibit superconductivity in the chan- 
nel. Most calculations of the proton critical temperature, 
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FIG. 3. A typical good fit to Cas A's rapid cooling for a 
1.4M star, built from the EOS of APR [!J with an envelope 
mass AA/n g ht = 5 x 10~ 13 M©. The two dotted curves, with 
indicated values of Tc, are to guide the eye. The three models 
have a proton 1 So gap from [25|] (the model "CCDK" in 
which results in the entire core being superconducting. The 
insert shows a comparison of our results with the five data 
points of [3] along with their la errors. 



T cp (p), are larger than T cn (p) at low densities. Proton densities than current models predict (see, e.g., Fig 9 in 



superconductivity suppresses the MU process in a large 
volume of the core at a very early age, reducing Lmu 
24J . In our analytical model, this reduction translates 
to a lower Lg and, hence, to a larger /. The analyti- 
cal model as well as our calculations reveal that proton 
superconductivity significantly accelerates cooling during 
transit and results in a large slope. This feature, essen- 
tial to account for Cas A's cooling rate, is illustrated in 
the right panel of Fig. [2] 

By varying the relevant physical ingredients, such as 
the density range of proton 1 Sq superconductivity, the 
shape of the T cn (p) curve, the chemical composition of 
the envelope, and the star's mass, many models can re- 
produce the average observed T e of Cas A. These models 
yield slopes ranging from ~ —0.1 (no rapid cooling and 
no constraint on Tc) up to —2. A typical good fit to 
the rapid cooling of Cas A is shown in Fig. 02 where the 
large slope results from the strong suppression of Lmu by 
extensive proton superconductivity. Fig. 0] demonstrates 
that the result Tc — 0.5 x 10 9 K does not depend on the 
star's mass, but that the slope during the transit is very 
sensitive to the extent of proton superconductivity. Mod- 
els successful in reproducing the observed slope require 
superconducting protons in the entire core. Although 
spectral fits [f| seem to indicate that Cas A has a larger 
than canonical mass (1.4 M©), a recent analysis |6( in- 
dicates compatibility, to within 3<r, with a smaller mass, 
1.25 M©. The need for extensive proton superconductiv- 
ity to reproduce the large observed slope favors moderate 
masses unless superconductivity extends to much higher 



14} for a large sample of current models) . 
The inferred T c ~ 0.5 x 10 9 K, either from Figs. [TJ 
O and S] or from Eq. (O, appears quite robust 
and stems from the small exponent in the relation 
Tc oc (Cg-Lg -1 ^ 1 ) 1 / 6 . Assuming L 9 is not very strongly 
affected by proton superconductivity, Tc is mostly de- 
termined by tc which, given the briefness of the rapid 
cooling transit phase, cannot be much smaller than the 
age of Cas A. On the other hand, the rapidity of the 
cooling, i.e., the slope s, is predominantly controlled by 
proton superconductivity reflected in the suppression of 
Lmu- We note, however, that if the proton gap ex- 
tends to much higher densities than current theoretical 
models indicate, suppression of Lmu could be larger than 
considered here. This suppression can be up to a factor 
~ 50; thus, Tc could almost reach 10 9 K [26| . 

Although we have assumed the minimal cooling 
paradigm |14| , our results remain valid in case some fast 
v process is allowed. Any of these processes is so efficient 
that the high observed T^° of Cas A implies that it has 
been strongly suppressed at a very young age (likely by 
pairing of one of the participating particles) and rendered 
inoperative 11] . In such a case, L„, prior to the onset 
of the neutron 3 P2 phase transition, was dominated by 
the MU and nucleon bremsstrahlung processes, exactly 
as in the minimal cooling paradigm, and the results of 
our present analysis are not altered. 

Our value of Tc deduced from Cas A's cooling is com- 
patible with the requirement Tc <^ 0.5 x 10 9 K inferred 
from the minimal cooling paradigm described in (l9| . 
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FIG. 4. Cooling curves with different masses and values of Tc 
as indicated. For the 1.9M Q star, AMi igh t = 5x 10 _11 M Q . 
For the other two masses shown, AMiight = 5 x 10 -13 M© . The 
assumed proton 1 So gap is the same as in Fig. [3] The slopes, 
at the current age of Cas A, are —1.4, —0.9, and —0.5 for the 
1.3, 1.6, and 1.9M© models, respectively: the decrease, with 
increasing mass, directly reflects the decrease of the core's 
fractional volumes in which protons are superconducting. 



We predict that Cas A will continue to cool for several 
decades with an almost constant rate, in contrast to two 
other interpretations Q; a relaxation from a recent sud- 
den energy release deep in the crust or a low level ac- 
cretion from a fossil disc that imply that the observed 
cooling is episodic and that the light curve should flatten 
out soon. Distinguishing among these possibilities and 
our model will be possible by continuous monitoring of 
the thermal evolution of Cas A. 
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